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The anomalous J/ip suppression in Pb-Pb collisions at 158A GeV observed recently by NA50 can 
be explained as due to the transition to a new phase of strong J/%j} absorption, which sets in when 
the local energy density exceeds about 3.4 GeV/fm^. 
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High-energy heavy-ion collisions have become the focus 
of intense research because of the possibility of producing 
a quark-gluon plasma during such collisions The 
suppression of J/ij} and "0' production has been suggested 
to probe the screening between a charm quark and a 
charm antiquark in the plasma |^J^. While J/^ and 
ip' suppression has been observed PFh| , the phenomenon 
can be explained by absorption models without assuming 
the occurrence of the plasma 

Recently NA50 observed that the J/ip and i/'' produc- 
tion is anomalously suppressed for Pb-Pb collisions at 
158A GeV The observation led to a flurry of ac- 

tivities. While the present author and others presented 
theoretical findings at the Quark Matter '96 Meeting in 
May, 1996 f 



other theoretical studies have since 
been put forth 1 2C J2^ . A central question is whether the 
anomalous suppression in Pb-Pb collisions arises from 
the absorption by comovers (produced hadrons), as pro- 
posed by or from the occurrence of a new phase 
of strong J absorbing matter (possibly a quark-gluon 
plasma), as suggested injjl6|-[l8| . The anomalous sup- 
pression of "0' in Pb-Pb 



L5[ and S-U collisions [[7|g] is an- 
other question which needs to be addressed in the context 
of any possible occurrence of the new phase. 

The study of these two questions requires the examina- 
tion of absorption by soft particles not in the new phase, 
which is beyond the scope of the schematic model of . 
Previously, a microscopic absorption model (MAM) was 
proposed which allows one to study the absorption by 
soft particles We shall use the MAM model to ex- 
amine the two central questions outlined above. We show 
that contrary to the conclusions of ]l9|-pT[, absorption 
by produced soft particles cannot explain the anomalous 
J /ip absorption in Pb-Pb collisions; that a class of mod- 
els with the feature of the occurrence of a new phase of 
strong absorption can explain the complete set of J/?/; 
data; and that 0' is already anomalously absorbed by 
produced soft particle not in the new phase, and the tran- 
sition to the new phase leads only to a small increase in 
tj}' absorption in Pb-Pb collisions. 

We envisage that a produced (quasi-bound) cc pair 
evolving into J /ijj and ip' states will collide with two 
types of particles. Collisions with baryons occur at high 
relative energies, and constitute the hard component of 
the absorption process. The absorption cross sections 



fabs('0-^) and (Jahsi'^P' N) at high energies are approxi- 
mately equal empirically |p2| , which can be understood 
by the Glauber picture of hadron-hadron collision [|TT| . 
The collisions of the cc with produced gluons or produced 
hadrons, which constitutes the soft component of absorp- 
tion, occur at low relative energies of about 200 MeV, the 
temperature of produced particles. The breakup thresh- 
old is about 640 MeV for J/-0 and only about 50 MeV 
for ip' . One expects the soft component to be small for 
J /ip but large for ip' . This is borne out by the data in 
Fig. 1 where Ba/AB is plotted (in logarithmic scale) 
as a function oi rj = A^/^[A -1)/A + B^/^{B - 1)/B, 
which is approximately proportional to the average path 
length passing through nuclei A and B [Q. The absorp- 
tion factor for the hard component is exp{— constantX77}. 
The absorption factor for the soft component is simi- 
larly exp{— constant' x rj} because the density of the pro- 
duced soft particles is proportional to the longitudinal 
path length passing through nuclei A and B (see pages 
374-377 of H]). p-A collisions involve the hard compo- 
nent while A-B collisions involve both the hard and soft 
components |pT[| . The magnitude of the soft component 
is indicated by a gap and a slope change between the A-B 
line and the p-K line in Fig. 1. A large soft component 
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for ip' absorption is indicated in Fig. lb by a large gap 
between the S-U point and the p-A line. A very small soft 
component for J /ip absorption is indicated in Fig. la. 



1 



as there is almost no gap and no slope change between 
the p-A line and the A-B line joining the 0-Cu, 0-U, 
and S-U points. The A-B line passing the 0-Cu, 0-U 
and S-U points is much above the Pb-Pb point in Fig. 
la, indicating that the soft component constrained to 
explain the data of 0-Cu, 0-S and S-U cannot explain 
the anomalous J/V' suppression in Pb-Pb collisions. A 
new source of absorption is suggested. 

We shall further confirm the above observations us- 
ing the microscopic absorption model (MAM) . Adopting 
a row-on-row picture in the center-of-mass system and 
following straight-line space-time trajectories of the cc's, 
the baryons, and the centers of the fireballs of produced 
soft particles, the differential cross section for J /-ip pro- 
duction in an A-B collision is 



da (b) 

a do 

.i/i, 



db. 



X 1- 



1 B 



l-r,(b-6Jaabs(^^) 



F{bA), (1) 



where TA{bA) is the thickness function of nucleus A, and 
the soft particle absorption factor F{bA) is 



N< 



F{bA)^J2 



a{n) 



iV>7V< 
1 > < j=i 



J2eM~0^ik^gtf^ + k^Ht^)} . (2) 



Here, Ny{bA) and A^<(b^) are the greater and the lesser 
of the (rounded-off) nucleon numbers AT/i(6^)CTi„ and 
BTsib — bA)<Jin in the row at bA with an NN inelas- 
tic cross section aim o,{n) = 2 for n ~ l,2,...,Af< — 1, 
and a{N^) — Ny — N^ + 1. The quantities tTabs(V'™)i 
the average relative velocity v^m and the average num- 
ber density pm per NN collision always come together 
so that we can use the rate constant fc^™ (m = g, h) 
to represent their product. The interaction time tf- (or 
t'j) is the time for a J /ip produced in collision j to over- 
lap with the center of the fireball of gluons (or hadrons) 
produced in collision i at the same spatial point. They 
can be determined from cc, 5, h production time tec, tg, 
th, and the freezeout time tf. The function 6 is zero if 
A = 1 or i? = 1, and is 1 otherwise. The expressions for 
■0' production can be obtained from Eqs. (|l|j^) above by 
changing -0 into ip' . 

We first study Ba/AB data without the Pb-Pb points 
and compare them with theoretical models, including the 
possibility of absorption by soft particles: (A) absorp- 
tion by baryons only, as i n (B) by baryons and pro- 
duced soft gluons, as in and (C) by baryons and 
produced soft hadrons (similar to ||lO(] but differing in 
details and methods of evaluation). Following [0, we 
use a^hsi^'N) = crabs (V'A^) and search for (Tabs(V'A^) as 



well as h 



and k^'g in Model B, 



and k^h and k 



Model C, fixing the time constants to have the plausible 



values tg = 0.1, th = 1.2, = 3, and tec = 0.06 (in units 
of fm/c). Different time constants will modify inversely 
the rate constants but will not affect greatly the product 
of the rate constants and their corresponding average in- 
teraction times. We take a Gaussian nuclear density for 
yl < 40 and a Woods-Saxon density for A > 40. The 
curves marked "without new phase" in Fig. 1 are calcu- 
lated with the parameters aahsitpN) = 6.94 mb in Model 
A, and (Tabs('0A^) = 6.36 mb in Model B and C. In addi- 
tion, k^g = 0.0956 c/fm and k^>g = 3 c/fm for Model B, 
and k^fi = 0.0493 c/fm fc^/^, = 3 c/fm for Model C. As 
kt(,g << k^ig in Model B and k^h << k^'h in Model C, 
we confirm that the soft absorption component is small 
for J/ip but large for 

When we include the Pb-Pb data point, no MAM cal- 
culations with a hard and a soft absorption component 
can simultaneously describe the whole set of J/tp data. 
This suggests that there is a transition to a new phase 
of strong J/^ absorption, when the local energy density 
exceeds a certain threshold. One can extend the MAM 
model to describe this transition. The energy density is 
approximately proportional to the number of collisions 
which has taken place at that point up to that time. We 
postulate that soft particles make a transition to a new 
phase of strong J/tp absorption if there have been Nc 
or more baryon-baryon collisions at that point at time 
tx- The quantity k^gtf^ -\- k^ht^j in Eq. (2) becomes 
^i>a^ij + f^i'htij + k^xtfj, where the new rate constant 
kt(,x describes the absorption of J/ip by the produced soft 
matter in the new phase. Here, tfj — tn + th ^ t^ where 
tn is the last NN collision time at that point and t^j 
is the time for a J/ip produced in collision j to over- 
lap with the center of the fireball of absorbing soft par- 
ticles produced in i in the form of the new phase, be- 
fore hadronization takes place. We vary N^ and fc^a^. 
Baryons passing through the spatial region of the new 
phase may also become deconfined and may alter their 
ijj-N absorption cross section. Accordingly, we also vary 
the effective absorption cross section, a^y^^{ipN), for ip- 
N interactions in the row in which there is a transition 
to the new phase, while crabs('0A^) remain unchanged in 
other rows. The curves marked "with new phase" in Fig. 
1 are obtained with the parameters N^ = 4, k^^ = 1 
c/fm, and a^^^^ipN) = 14 mb. The location where the 
slopes of the curves changes sharply is sensitive to Nc- 
The absorption of J/ip saturates for k^x ^ 1 c/fm for 
which the absorption is still slightly insufficient to ac- 
count for the total J/ tp absorption; additional absorption 
with o'fbs(V'-^) = 14 mb is needed to give Pb-Pb result 
to agree with experiment. The results for the cases with 
the new phase agree with the whole set of J/ ip data. We 
note that k^x » k^g, k^h- 

The above MAM results can be expanded to obtain 
J/ip and ip' yields as a function of impact parameter b. 
The theoretical ratio Bajf^/a"^^ {DY)2.9^4,.5 for Drell- 
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Yan cross section in the interval 2.9 < M^+^- < 4.5 GeV 
can then be expressed in terms of Baj^^^/aP^{DY)2.9-4.5 
for pp colhsions, which was determined by NA51 to be 
44±3 [|l5|. NA38 and NA50 used a geometrical model 
to obtain a relation between Et and b which can 

be used to transform the MAM results from a function 
of 6 to a function of Et- In Fig. 2a we show theoretical 
Bajf^/a^^{DY)2.9-4.5 as a function ofE^ for S-U which 
agree with experiment within 5-10%, whether we assume 
a new phase or not. On the other hand, good agreement 
with Pb-Pb data is obtained only when we allow for the 
transition to a new phase of strong J/ip absorption (Fig. 
2b). 
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Fig.2. Ba^^^l a[DY)-2.9-A.b as a function of E^ for 
(a) S-U collisions, and (6) Pb-Pb collisions. Data are 
from NA38 [| and NA50 

The results in Fig. 16 indicate that ij]' is already anoma- 
lously absorbed by soft particles in S-U and Pb-Pb col- 
lisions, even without the new phase. The occurrence of 
the new phase with k^ix = 3 c/fm leads only to a small 
additional ij)' absorption in Pb-Pb collisions (Fig. \b). 
We can calculate B'a{ip')/Ba{J/tp) and compare with 
experiment. Fig. 3& shows that the experimental data 
for Pb-Pb collisions is consistent with the assumption of 
transition to the new phase. For S-U collisions, the theo- 
retical results agree with experimental data for small and 
moderate values of Et, but deviate from experimental 
data for large Et (Fig. 3b). As large transverse ener- 
gies involve greater weights for the major axis of the U 
nucleus to lie along the beam direction, the deviation of 
B'a{ip')/Ba{J/4>) in S-U collisions at high Et may be a 
deformation effect which can be tested experimentally by 
studying S-Pb collisions. The deformation effect can be 



utilized to study matter in the new phase by focusing on 
high Et events in U-U collisions. 
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Fig.3. B'a{il^')/Ba{'4}) as a function of Et for (a) 
S-U collisions and (b) Pb-Pb collisions. The data for 
p-W and p-U collisions are also included. Data are 
from NA38 lilTtl and NA50 [pi. 



What is the threshold energy density for the transi- 
tion? We consider a row with a transverse area (7 in 
and with collision points separated longitudinally by d/j, 
where d = 2.46 fm is the internucleon spacing and 
7 = y/s/2mnuc\ = 9.2 is the Lorentz contraction factor at 
the Pb-Pb collision energy. Consider also Nc number of 
NN collisions at each collision point. Each NN collision 
leads to a fireball of soft particles with dn^^ /dy ^ 1.9 at 
the Pb-Pb collision energy |^^. At the time t — d/^ af- 
ter these Nc collisions, the produced soft particles leaving 
the fireballs at one collision point to the adjacent collision 
point will be compensated by soft particles arriving from 
the fireballs of the adjacent collision point, and a steady- 
state initial energy density is reached after t = d/7 
(before longitudinal expansion), with the energy den- 
sity given approximately by Nc{dn^-^ /dy)mt/{aind/^), 
where mt = 0.35 GeV is the pion transverse mass. Hence, 
for Nc — A from the present study, the threshold en- 
ergy density for the new phase is Ec ~ 3.4 GeV/fm"^, 
which is close to the quark-gluon plasma energy density, 
Ec ^ 4.2 GeV/fm'^, calculated from the lattice gauge re- 
sult Cc/T^ - 20 [|6| with Tc - 0.2 GeV. Therefore, it is 
interesting to speculate whether the new phase of strong 
J/f/; absorption may be the quark-gluon plasma. In the 
equilibrated or non-equilibrated quark-gluon plasm, J/^p 
production is expected to be greatly suppressed P,p7[. 
Furthermore, when baryons pass through the region of 
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the new phase, the baryon matter may be deconfined. 
The total cross section between a cc system and a baryon 
system is substantially enhanced when the quarks in the 
baryon are deconfined psf . 

Our results for the cases with soft particle absorption 
differ qualitatively from those of ]l9|-pl[|. We would like 
to mention some of the reasons for the differences. In 
Refs. | p^ , pO| , the soft particle absorption factor is ap- 
proximately the form exp{— constant x Exib)}. This is 
obtained in by assuming the density of produced soft 
particles to be n{b) = constant x ET{b). Ref. as- 
sumes n{b) = constant x G{b)ET{b) and finds G{b) to 
be independent of b. However, Exib) contains contribu- 
tions coming from the whole volume of NN collisions of 
the participant nucleons, which decreases significantly as 
b increases. When this volume V{b) (in the CM. sys- 
tem) is taken into account, the proper relation should be 
n{b) = constant x ET{b)/V{b). The additional V{b)-type 
dependence in n{b) will modify the results of [ p^j20| . Fig. 
2 of the numerical cascade model of gives a resultant 
J/ absorption factor for J/i/j absorption due to baryons 
in S-W collisions to be smaller than that in p-W colli- 
sion at the same impact parameter of 2 fm. The authors 
of should check whether these results are consistent 
with intuitive understanding and analytic results for ab- 
sorption due to baryons. Clearly, much work remains to 
be done to resolve the differences. 
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